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Abstract 

The renormalised value of {4> 2 ) is calculated for a massless, confor- 
mally coupled scalar field in the Hartle-Hawking vacuum state. This 
calculation is a first step towards the calculation of the gravitational 
back reaction of the field in a black cosmic string spacetime which 
is asymptotically anti-DeSitter and possesses a non constant dilaton 
field. It is found that the field is divergence free throughout the space- 
time and attains its maximum value near the horizon. 



1 Introduction 

A useful question to ask when one studies quantum fields in General 
Relativity is the following: Given that all matter is inherently quantum in 
nature, will quantum effects remove the singularity at the centre of a black 
hole spacetime? To answer this question using a scalar field and semi-classical 
perturbation theory, one must first calculate the expectation value (T^) 
where is the stress-energy tensor operator of the scalar field <ft. In this 
paper (0 2 ) is computed in preparation for the computation of (T^) which is 
used as the source term for the Einstein field equations: 
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R P v - )^Rg^ - ^9tiu = 8vr(T^). (1) 

A review of handling quantum fields in the presence of strong gravitational 
fields can been found in articles by Wipf jjj] and DeWitt |2| and in books 
by Birrel and Davies and Wald (|. 

((f) 2 ) for massless fields has been computed for both the interior and 
exterior of a Schwarzschild black hole H] ||. These calculations have also 
been extended by Anderson to accommodate massive fields in general spher- 
ically symmetric, asymptotically flat spacetimes j7|. A method has also been 
developed by Anderson, Hiscock and Samuel || |{J to calculate the expecta- 
tion value of the stress-energy operator, (T^ u ) in spherically symmetric, static 
spacetimes. They use this method to calculate (T M „) in Schwarzschild and 
Reissner-Nordstrom geometries. (T^) has also previously been computed by 



Howard and Candelas |10[ in Schwarzschild spacetime. 

The Kerr spacetime has also been studied with fields propagating 
in this geometry. Frolov [ll] has calculated (0 2 ) for massless fields on the 
event horizon pole of a Kerr-Newmann black hole as well as deriving an 
approximate expression for (T^) near the horizon with Thorne Massive 



fields in the exterior geometry have been studied by Frolov and Zel'nikov JT3 

Various works on back reaction effects of quantum fields on black hole 
geometries have been produced. For Schwarzschild and Reissner-Nordstrom 
spacetimes this includes the work of Hiscock and Weems ||14]| , Bardeen and 



Balbinot [|1J |T(|, and York [|17| who used Page's analytic approxima- 



tion Jl8| for (T^) in Einstein spacetimes for conformally invariant fields. 
More recently Hiscock and Larson |19| have extended their analysis to the 
Schwarzschild interior and calculated back reaction effects on curvature in- 
variants. 

Most work in this field has been done in the context of spherical or 
oblate symmetry. We wish to extend the analysis to other symmetries and 
ultimately ask if the above results are general or are specific to the particular 
symmetry chosen. For example, does the presence of the quantum field have 
the same effect on the curvature growth and on the anisotropy of a black 
hole interior for all symmetries? If curvature invariants are weakened for all 
systems studied then one can say with some confidence that quantum effects 
may remove the singularity. 
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1.1 Background geometry 

The system studied here will be that of a massless Klein-Gordon field 
with conformal curvature coupling propagating in the spacetime generated by 
a straight black string. This background is chosen for several reasons. First, 
it possesses cylindrical (as opposed to spherical) symmetry and secondly, in 
the context of cosmic strings it represents a system which may physically 
exist in the universe. It could also be argued that sufficiently close to a 
black string loop the spacetime will possess this type of geometry. There 
has also been a revived interest in anti-DeSitter spacetimes in the context 
of conformal field theories. The metric chosen is that developed by Kaloper 



pOf and Lemos and Zanchin uB which admits a line element of the following 



form: 



= - w -=i+$4'-gi- pay dt&p 



+ 



"I 2 2 2(3Q-A/) (30-Af)4Q2- + °^ ' P 2 dz 2 (2) 

01 P op + (n + Af)c*2 p 2 

Where M, Q, and J are the mass, charge, and angular momentum per unit 
length of the string respectively. Q is given by: 



The constant a is defined as follows: 

« 2 = -^A (4) 

with A, the cosmological constant, negative giving the spacetime its asymp- 
totically anti-DeSitter behaviour. The spacetime has a well defined time-like 
killing vector field with respect to which modes can be defined. 

Much interesting work has been done regarding effects of quantum 



fields in the 3D BTZ |2[ black hole H,&BH,|g,[g|]. The system 
here however does not dimensionally reduce to the 3D BTZ black hole whose 
four dimensional counterpart (without charge or angular momentum) is given 
for reference as f29fl 

dv 2 

ds 2 BTZi = -(a 2 r 2 - 8M)dt 2 + - + r 2 d V 2 + dz 2 . (5) 
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This is due to the fact that the dilaton field is non zero and non constant in 
the corresponding three dimensional action. 

In the spacetime considered here, both charge and angular momen- 
tum are zero yielding the following for (0) 



ds 2 = -{a 2 p 2 - — ) dt 2 + - 2 d 2 P _ 4Mx + p 2 V + a 2 p 2 dz 2 . (6) 



ap [a*p< - ^) 



It can immediately be seen from (|J) that the spacetime behaves as anti- 
DeSitter as p —>■ oo and therefore the spacetime is not globally hyperbolic. 
This fact requires boundary conditions to be imposed to control the flow of 
information in and out of the time-like infinity. It has been shown |3(| that 



three natural boundary conditions arise in ADS4 spacetime and the condition 
used here is that of a "transparent" boundary condition. 



2 Green Function and (0 2 ) 

In this section a Euclidean space approach is used to calculate the Green 
function. The calculation is similar to that done by Anderson who cal- 
culated (0 2 ) in a Reissner-Nordstrom spacetime however, the method is ex- 
tended here to a system which is neither spherically symmetric nor asymp- 
totically flat. 

The Euclidean space method amounts to making the transformation 
t — > it in metric (@). The quantity ((f) 2 ) is then defined as 

(</> 2 > = limC^x'), (7) 
where Ge(x, x') is the Euclidean space Green function satisfying the equation 

[□ a - m 2 - ZR(x)} G E (x, x') = (8) 

The mass term in the Klein-Gordon operator will be set to zero in order to 
comply with the transparent boundary conditions mentioned earlier. 

The presence of three killing fields allows the right-hand-side of the 
above equation to be expanded in terms of cylindrical functions as follows: 
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with 



8{t-t>) = T « 

n=— oo 
1 oo 

- y 



= in27rT(r-r') 



% - <p') 

- z') 



(10) 



Z=— oo 



— / dke tk{z - z '\ 
2tt 



Where T is the temperature of the field. The Green function can similarly 
be expanded giving 

rp OO OO /»00 

G £ ( X;X ') = _i_ e m2 ^ T - T>) \ dke^ z -^ X {p,p')- (H) 



Z=— oo 



The function x(PiP') must satisfy the following equation 



d 2 x(p,p') 



dp 2 



2 „2 



4M 



a p — 



+ 



c?x(p, P') 



-x(p,p r ) 



,2 /-»,2 /-i2 



+ 



ap y dp 
n 2 4vr 2 T 2 



4 oTp - 



M 



a 2 p 2 - 



p^ crp 

A solution is assumed of the form 

X(p,p0 = C^ 1 (p < )^ 2 (p > ), 



ap 
p 2 a 



;i2) 



(13) 



where p< and p> represent the lesser and greater of p and p' respectively. In- 
tegrating across the 5 function gives the Wronskian normalisation condition: 



C 



dp 



dp 



ap 2 (a 2 p 2 - 



(14) 



The asymptotic behaviour of the solution can be found by studying 



the solutions to ( |T2"D in the appropriate regimes. For p = oo the solution has 
the form: 



Whereas near the horizon yp 

^2 



(4M) 



1/3 



solutions are found to behave as: 



ap dp 
o 3 p 3 -4M 



(16) 



Where k = 2tcT. It can easily be seen that diverges at infinity and ^2 is 
divergent at the horizon. The general solution with the correct asymptotic 
behaviour can be found by ansatz. Such a solution takes the form 



#1 
^2 



1 



f p2 a 3 X dp, 
g J r r V 900333 



1 - r p2 a 3 Xd jZEi 



X is a function of p which evolves according to the equation: 



(17) 



X 2 



a 2 p 2 



AM 
a 3 p 3 



l 2 + 



k 2 



Q 



+ 



K 2 n 2 



a' 



1 - 



AM 
ot A p A I 



+ T - 



36M 2 

+ ^7 



a 6 /? 6 



+ 



~X 



a" 



P 2 + 





4M \ 


0- 


a 3 p 3 J 


6M 


12M 2 ' 


a 3 p 


a e p 4 _ 



1 2M 
2 P+ 



16M : 



21 



(18) 



+ 



AM 8M 2 



a 3 p 



a 6 p 4 



Which can be obtained by substituting ( |TTD into (§). Substituting ( |T7| ) into 
fll4f ) gives C = I for all mode functions. 



A point splitting algorithm developed by Christensen |H| will be 
used to renormalise the field. In this technique, one chooses the points x and 
x' to be nearby points in the spacetime before the full coincidence limit is 
taken. It is convenient to have the points take on equal values of p, (p and z 
so that the coordinate separation is given by e = r — r'. The unrenormalised 
Green function now takes on the form: 



rp 00 00 p 



cMi(p)^ 2 (p). 



(19) 



It should be noted at this point that there exists a superficial ultra-violet 
divergence over I and k in the above expression. This divergence can most 
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easily be eliminated using a similar technique as Candelas || and Anderson 
0. It is noted that as long as r ^ t' any multiple of 5(t — r') can be added 
to ([[£]). Substituting X from ( |I8D in the large Z and limit and subtracting 



this term from (|T^) the logarithmic divergences can be eliminated giving the 



following expression: 



Ge(x,t,t>) = — e™ [ J2 J_JkMp)Mp) 

\ 



n=—oo \l=— oo ' 



(20) 



2.1 Renormalisation 

To calculate the renormalised value of (0 2 ) a point splitting technique 
will be used. The DeWitt generalisation to Schwinger's expansion is used as 
an approximation for the Green function. This term will then be subtracted 
from (|20|) and the x — > x' limit will be taken along the shortest geodesic 
separating the points. The DeWitt-Schwinger counter-term is given by: 

1 1 a'^a'< v 

G(x,x') = t-^ + ^R,u (21) 

8n z (T 9bn z a 



where a is the "world function" of Synge |32] which is equal to half of the 
square of the geodesic distance between two points. The points in this case 
will be t and r'. 

It can be shown, by geodesic expansion (see appendix A), that the 
world function in the spacetime considered here takes on the form 



2 V up J 24 a 3 p 5 



4 i n(J>\ 



1 (a 3 p 3 - 2M) (a 3 p 3 - AM) 
6 a 3 p 5 



^ = v-r ^ e 3 + (e 5 ), (22) 



(q 3 p 3 + 2M) (q 3 p 3 - 4M) 2 

= a"! e + <A e 
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Using these expressions the DeWitt-Schwinger counter term is equal to 

1 1 (a 3 p 3 + 2M) 2 

4tt 2 (a 2 p 2 - —\ e 2 4871 " 2 a <° 3 ( a3 '° 3 ~~ 

(23) 



48vr 



The first term in (^) can be rewritten in a more convenient way by 
using the Plana sum formula [p3] as used by Anderson. The formula is 



1,, v, f°,, x, , [°°dt [f{n + it)-f{n -it)\ 



v i poo n 

J2f(n) = -f(n )+ / /( S )efa + l / 

n=nn ^ ^"O ^0 



e 27Tt _ 1 



(24) 



With this, the first expression in the counter term can be written as 

1 



4vr 2 (a 2 p 2 - e 2 



OO o 



Ecos(n/te) n k -. r-, (25) 



An 2 [a 2 p 2 -^)^ " ' 48vr 2 (a 2 p 2 - ^) 

which gives, for the entire counter-term in the e — > limit: 



r K y 

^counter / \ / , 



OO n 

n k 



47T 2 (« 2 p 2 - ^) ^ 487r 2 (aV-^ 
1 (a 3 p 3 + 2M) 2 3a 2 



48vr 2 ap 3 (a 3 p 3 - AM) 48vr 2 



(26) 



It can be shown that the second and third terms in (p6|) , which normally 
diverge at the horizon, will cancel each other out on the horizon when T is 
equal to the black hole temperature. That is, when 

T = £~(4M)*. (27) 
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Therefore, the renormalised expression for (0 2 ) in the Hartle-Hawking vac- 
uum state is 



T 



\ 



8vr 2 ( Yl 



-oo l=—oo 



+ 



4tt 2 / 



E 



n k 



dk^ 1 (p)^ 2 (p) 

i 

f 2 



/2 + S 

3a 2 



a 2 p 2 - *M 
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n=l 



48vr 



2 ' 



48tt 2 / 192tt 2 / 

Where, for convenience, the following notation has been used 

4Af 



ap 
(28) 



f = 9oo 

r - 



2 2 

a p 



ap 



df 
dp 



(29) 



3 Calculation of (0 2 ) 



In this section the value of (<j) 2 ) will be calculated from ( 
solution is found by iteratively solving for the function X using 
mode functions with the lowest order term defined as: 



h 



a 2 p 2 



l 2 + 



+ 



K 2 n 2 
a 2 h 



+ 



9 



36M 2 £R , 
R a H irh, 



and h given by 



h 



AM 
a 3 p 3 



. The 
in the 



(30) 



(31) 



With this choice, the lowest order term will be valid at both the horizon and 



as p approaches infinity. This can also be verified by comparison with ( 15 ) 
and flTEl). 

It is most convenient to do the integration over k first followed by 
the sum over / and finally, the sum over n. This scheme leads to analytic 
solutions to both the integrals and the / sums. 

In the solution, the sums and integrals have the following form: 



EE 



dk 



P 2 + S + S£ + VM) 



p ■ 

2 



(32) 
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where Vq(p) is a function of p only. Such integrals are known and are given 
by 

dk „ a , . 

C- - (33) 



(^ 2 + S + W + ^o(p)) 5 (/ 2 + ^ + ^o(p))' 



2 



where C is a fractional constant which depends on the particular value of 
the integer p. The resulting sums over / can now be done analytically by the 
standard contour integration: 



1 NT^ / „ . , , 7T COt 7T I 

Residues of 



2 



(34) 



This is valid as the denominator never becomes singular at integer values of 
In the expansion of (0 2 ) there appears a term of the form 



which arises when the lowest order term is combined with the counterterm in 
(|20| ) and integrated over fc. This is the only term which is singular at I — 0. 
The / = term is therefore removed from the sum and redefined as to remove 

a spurious divergence of the form In (^ 7! '\J^i^- + K)(pTj which is present in 

the same term once the sum from / = 1 to oo has been calculated (how this 
term arises is shown in more detail in appendix B.) This redefinition can be 
performed since the counterterm introduced in ( ^0|) may be any quantity 
and the particular value shown there is exact only in the large I limit. 

The final sum over n can now numerically be shown to converge by 
computing the values of (<f> 2 ) for large n. It can also be shown that the n = 
mode makes no contribution to the sum. 

The boundary values of the modes must be evaluated before numer- 
ical integration of the mode equation can be done. The value of the mode 
functions at infinity can easily be seen by studying the mode equation in the 
asymptotic region and therefore only the value at the horizon is left to be 
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determined. At the horizon, there are many quantities in the expansion of 
((f) 2 ) which are inversely proportional to some power of the metric function 
/. By performing an expansion of ((f) 2 ) in the quantity 5 = p — pn, where pu 
is the horizon value of p, one can show (although the procedure is lengthy) 
that all terms with 5 raised to some power in the denominator cancel at the 
horizon. In appendix B it is also demonstrated how terms which normally 
make a dominant contribution to ((f) 2 ) cancel here. It should be noted that 
the horizon value is directly proportional to the value of the cosmological 
constant. This is due to the fact that the spacetime is an Einstein spacetime 
with = 3a 2 g^. 

((f) 2 ) was computed for the conformally coupled, massless case and 
the result is shown in Fig.l. It can be seen the the maximum value of ((f) 2 ) 
occurs near, but not at, the horizon. This behaviour is analogous to the ex- 
treme Reissner- Nordstrom case J7| . This is because, as shown in the appendix 
and earlier, most contributions to ((f) 2 ) at the horizon vanish. However, near 
(but not on) the horizon terms with a 1/ f behaviour make a large contribu- 
tion. For large p most terms in the field expansion vanish and therefore ((f) 2 ) 
approaches a value which is dominated by the last two terms in (p8| ) . 




1.59 1.79 
Radius 

Figure 1: ((f) 2 ) for the cylindrical black hole spacetime. The value of ((f) 2 ) 
has a small but nonzero value at the horizon and attains a maximum away 
from the horizon. 
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5 Appendix A 

In this appendix it will be demonstrated how to calculate the world 
function, a, using the method of geodesic expansion. Equations here involve 
quantities which are defined at different spacetime points. For a brief review 
of handling bitensors the reader is referred to Christensen [|3l| and Synge 



0- 



Let Pi(x) and P2(x') be two points in the spacetime close enough together 
such that they are connected by a unique geodesic. The geodesic equation 

+ C-7T-7T = (36) 



ds 2 a(3 ds ds 
yields the power series 

*' + U»' ds - hi, pl U a U^ds 2 + 1 (lY^T^UWU 5 ) ds 3 



X M « x^ 



where 



and 



^^IPlPds 3 + ... (37) 



U» = ^ (38) 
ds 



U^ds = dx" w e = (r - t')5%. (39) 



The last term in equation (|3"T| ) is zero since the spacetime is static. ( P7[ ) can 
be inverted and used in the definition [B31 



2a{x, x') = ds 2 g^'U^'U"' (40) 
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giving 



+ ^eW^ + jr^rrr^r^ Vl/ . (4i) 

Calculating the Christoffel symbols using (|6|) and noting ([39]) it can be shown 
that this expression reduces to the one in (0). This expression has the 
same functional form as that of the world function for a static, spherically 
symmetric spacetime if the coordinate separation there is also chosen as 
fl39|). This is due to the fact that the £ M vectors eliminate any dependence of 
a(x,x') on #22 and g 33 . 



6 Appendix B 

In this appendix the dominant terms of (<ft 2 ) at the horizon will be 
calculated. This is useful since the value at the horizon needs to be evaluated 
as a starting point for the calculation of the mode functions. This calculation 
is also useful as it provides insight as to how the n counter-term acts to 
regularise the field. 

At the horizon, the dominant terms in the field expansion are given 

by 




(42) 

For the moment, we choose to ignore the I = term and concentrate on the 
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first expression in (^) which can trivially be re-written as 

/ \ 



T 

2r? 



oo oo 



n=— oo (=1 " v 



In 



W 1 + J ^ — / 



47T 



T 



p 2 a\fh 



(43) 



The product in the above expression is well known yielding the following for 

m 



rp OO 

— V 

Air 2 



4tt 2 ^ p 2 aVh 



n=— 00 

In I 2vr 



In 2 sinh 



K 2 n 2 
a 2 h 



K 2 n 2 



a 2 h 



+ Vo(p) 



(44) 



For very large n or (as in the case here) very small h this becomes 



T 



At: 2 



00 1 

E 



/t|n|7r / 2fi;|n|7r 
— m 



a 



a 



If we define the I = term to cancel out the second term in 
in the text, the resultant expression gives 



2vrT 
'4^7 



E 

n=0 



(45) 

as described 
(46) 



which is cancelled by the n sum in the counter-term. This leaves a small 
constant contribution to (<p 2 ) at the horizon. 
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